On the critical exponent of a quantum noise driven phase transition: the open system 

Dicke-model 
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The quantum phase transition of the Dicke-model has been observed recently in a system formed 
by motional excitations of a laser-driven Bose-Einstein condensate coupled to an optical cavity 
[T]. The cavity-based system is intrinsically open: photons can leak out of the cavity where they 
are detected. Even at zero temperature, the continuous weak measurement of the photon number 
leads to an irreversible dynamics towards a steady-state which exhibits a dynamical quantum phase 
transition. However, whereas the critical point and the mean field is only slightly modified with 
respect to the phase transition in the ground state, the entanglement and the critical exponents of 
the singular quantum correlations are significantly different in the two cases. 

PACS numbers: 37.30.+i,05.30.Rt,42.50.Nn 
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I. INTRODUCTION 

Experiments with ultracold atomic gases in optical 
fields laid down a new path to discover strongly corre- 
lated many-body quantum systems. In particular, the 
high degree of control over the interaction parameters al- 
lows for using atomic systems as quantum simulators of 
generic theoretical models 0. Central to these efforts 
lies the possibility of observing quantum phase transi- 
tions (QPT). At effectively zero temperature (T = 0), 
by tuning an external field acting on the system, it can 
be scanned through a quantum critical point which sep- 
arates regions with different symmetries in the ground 
state. One celebrated example is the QPT from a su- 
perfluid to a Mott insulator in the Bose-Hubbard model 
[S, that was realized with a gas of ultracold atoms in an 
optical lattice f^. Additional quantum phases appear in 
this system when dipole-dipole interaction is present [5]. 

A fundamental question is how quantum phase transi- 
tions are infiuenced by non-equilibrium conditions. The 
ordinary way to prepare a stationary system out of equi- 
librium at T = can be illustrated by a BEC in a rotat- 
ing trap. It undergoes the vortex formation QPT above a 
critical angular velocity [6j . External driving can impose 
that only a certain subset of states in the Hilbert space, 
those having a given moment of inertia in the previous 
example, be populated. Similar effect has been described 
for a spin chain in ring geometry: it can manifest crit- 
icality while being confined into the subspace of energy 
current carrying states |7j . In both examples the system 
is effectively Hamiltonian. 

One can go beyond the effectively Hamiltonian sys- 
tems by adding external non-equilibrium noise on critical 
states. It was shown that the 1// noise, ubiqitous in elec- 
tronic circuits, preserves the quantum phase transition in 
the steady state of a system, moreover, it gives a knob 
to tune the critical exponent by the noise strength [5]. 
This is in sharp contrast with the well-known effect of 
thermal fluctuations that destroy quantum critical corre- 
lations. In a more general level, reservoir engineering is 
a route towards designing specific noise sources in a dis- 



sipation process which leads to pure many-body states in 
the dynamical steady state. An example is a lattice gas 
immersed in a BEC of another species of atoms [9^ , which 
serves as a zero-temperature reservoir of Bogoliubov ex- 
citations. The resulting dissipative Bose-Hubbard model 
exhibits a dynamical phase transition between a pure su- 
perfiuid state and a thermal-like mixed state as the on- 
site interaction is increased [10]. Note that this method 
for the preparation of strongly correlated quantum states 
makes dissipation to be a resource for quantum simula- 
tion |11| and universal quantum computation |12j . 

In this paper we will consider the bare electromag- 
netic vacuum at T = as a reservoir and its effect on a 
Dicke-type Hamiltonian system which is known to pro- 
duce a singularity of the ground state [13] . Placed into a 
dissipative environment, the system evolves irreversibly 
into a steady state which is a dynamical equilibrium be- 
tween driving and damping. The intrinsic noise accom- 
panying the dissipation process is in accordance with the 
dissipation-fiuctuation theorem. Even in this very nat- 
ural case of non-equilibrium, the loss does not destroy 
quantum criticality. But what is the relation of the criti- 
cality expected in the steady-state to that of the ground 
state in the closed Hamiltonian system? 
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FIG. 1: Self-organization phase transition of a BEC in a cav- 
ity. Below a threshold in the transverse driving field (left) the 
condensate is quasi-homogeneous, and there are no photons 
inside the cavity. Above threshold (right), a standing matter 
wave of period A appears that scatter photons into the cavity. 



Our specific example is the self-organization phase 
transition of laser-driven atoms in an optical resonator 
[TH - fro] . The laser impinges on the atoms from a direction 
perpendicular to the resonator axis (see Fig.fTl). Below a 
critical value of the pump intensity, the spatial distribu- 
tion of the atoms is homogeneous along the axis and the 
mean cavity photon number is zero, since the photons 
scattered by the atoms into the cavity interfere destruc- 
tively. Above a threshold pump power, there appears a 
wavelength-periodic modulation of the distribution, from 
which laser photons can be Bragg-scattered into the res- 
onator. Spontaneous symmetry breaking takes place be- 
tween two possible solutions for the cavity field phase and 
the atomic distribution. This is a non-equilibrium phase 
transition, which has an experimentally accessible T = 
limit if the atomic cloud is represented by a Bose- Einstein 
condensate. The phase diagram has been experimentally 
mapped by Baumann et al. [1]. 



II. OPEN SYSTEM DESCRIPTION 

Consider the dispersive coupling of a one-dimensional 
matter wave field '^{x) to a single cavity mode a in the 
transverse pump geometry shown in Fig. [I] The disper- 
sive limit appears when the laser pump is far detuned 
from the atomic resonance (A^ = lu — uja exceeds the 
atomic linewidth 7 by orders of magnitude) . In the frame 
rotating at the pump frequency w, the many-particle 
Hamiltonian reads 



2 m dx^ 



H/h = -Aca^a+ f *^(a;) 
Jo 

+ Uq a' a cos^ (kx) + irjt cos kx{a^ — a) 



^{x)dx. (1) 



The detuning Ac = w — wc defines the effective photon 
energy in the cavity. Atom-atom s-wave collisions are ne- 
glected, the length of the condensate along the cavity axis 
is L. The atom-light interaction originates from coherent 
photon scattering. The absorption of a cavity photon and 
stimulated emission back into the cavity gives rise to the 
term proportional to Uq = g^ / /S.a- The coherent redis- 
tribution of photons between the pumping laser and the 
cavity mode results in an effective pump with amplitude 
r]t = Vtg/ /S.A. Note that this term describes the external 
driving of the system, and the explicit time-dependence 
in the optical frequency range, due to the laser field, has 
been eliminated by the transformation into the rotating 
frame. The remaining frequencies are in the kHz range 
of the recoil frequency ujj^ = hk'^ /2m. 

The critical behaviour can be described in a subspace 
spanned by two motional modes, i.e.. 



^{x) 



Tl' 



c\ cos kx . 



(2) 



the Hamiltonian of the system formally reduces to that 
of the Dicke model [20] . Originally, it was introduced to 
describe the dipole coupling of N two level atoms to a 
single quantized field mode [21J. It is known for a long 
time that the Dicke model can exhibit a thermodynamic 
phase transition at finite temperature [22], and a quan- 
tum phase transition at zero temperature |13| between 
an unexcited normal phase and a superradiant phase, 
where both the atoms and the mode are macroscopi- 
cally excited. There is a renewed interest in studying 
the zero-temperature properties of this system with par- 
ticular respect to critical entanglement [23| . finite-size 
scaling [231 US] and quantum chaos [2S]. Dimer et al. 
proposed a realization of the Dicke model with photon 
loss by means of multilevel atoms coupled to a ring cavity 
mode via Raman transitions 27J. Another collective spin 
model exhibiting dynamical QPT, the Lipkin-Meshkov- 
Glick model was constructed in cavity QED systems '28] . 
We consider a single dissipation channel which is the 
photon leakage through one of the mirrors. The cor- 
responding dissipation process can be modeled by a 
Heisenberg-Langevin equation for the field amplitude a, 
which includes a loss term with rate k and a Gaussian 
noise operator ^(t), 



dt 



-i\a ,li\ — Kfl -|- ^ . 



(3) 



with the bosonic annihilation operators cq and c\. With 
the closed subspace constraint CqCq -V c[ci = N imposed. 



The effect of continuous weak measurement of the photon 
number is described by the same equation. The noise op- 
erator ^ has zero mean and its only non- vanishing correla- 
tion is {^{t)^^t')) = 2n5{t - t') at T = 0. For finite tem- 
perature, other correlations would appear proportional 
with the thermal photon number. The given second- 
order correlation expresses the fluctuation-dissipation 
theorem. The noise operator can be seen as a necessary 
source for maintaining the commutation relation and gen- 
eral algebraic properties during the time evolution. 

This equation applies to the decay of a single uncou- 
pled harmonic oscillator. However, in the present case, 
the field mode interacts with the matter wave field. The 
decay process takes place in the optical frequency range 
which is many orders of magnitude above the charac- 
teristic frequencies of the interaction (comparable to the 
recoil frequency in the kHz range). Therefore, the cou- 
pling to the other atoms has a only a negligible influence 
on the decay of the photon mode and Eq. pi holds for 
the interacting system |29j . But one must keep in mind 
that the original problem is intrinsically time-dependent, 
and thus there is an energy current from the laser into 
the reservoir through the system. 



III. STEADY STATE VS. GROUND STATE 

Taking the N ^ 00 limit, the mean- field approach, 
which can be adopted both in the lossy and lossless cases 
of the Dicke- model (see Appendix A), gives an adequate 
approximation of the true steady-state and ground state. 




where < /3o < 1. The mean- field amplitudes |aoP, /Sg 
are plotted by thick lines (right scale) against the pump- 
ing strength y = ^/2N rjt normalized to the critical value 
j/c- The critical point separates two phases with differ- 
ent symmetries. In the normal phase ao = /3o = 0, in 
the superradiant phase there are two mean-field solutions 
corresponding to ±ao and 71 = ±/3o- The critical point 
He = depends on k. However, the mean-field solutions 
ao and /3 for the steady-state and for the ground state, 
expressed as a function of y/yc, overlap. The perfect 
overlap is not a general property, but there is always a 
smooth connection so that they coincide for k —>■ 0. This 
is because the mean-field solution is dominated by the 
identical Hamiltonian part of the dynamics. In contrast, 
the correlation functions of quantum fluctuations may 
signify different kind of second-order phase transitions, 
because the quantum noise source in the open system 
is dissipation (or the back-action of the measurement), 
which is completely missing from the closed system. 

For the matter field we consider the excitation mode 
Sb, which is orthogonal to the condensate wave function. 
In the normal phase, db = Sci, while in the superradiant 
phase Sb — —PqSco + \/l — PqSci. The figure shows the 
number of incoherent photons (da'Sa) (left scale, thin 
solid red line) and the condensate depletion (Sb^Sb) (thin 
solid blue line). The same physical quantities are plot- 
ted with thin dashed lines. There is a divergence for both 
cases, however, the exponent is —0.5 for the ground state 
(usual mean- field exponent) whereas it is —1 for the the 
steady state (see log- log scale in Fig. Isl) . This difference 



FIG. 2: Comparison of the criticality in the ground state of 
the closed, and in the steady-state of the open Dicke models. 
The mean fields and the incoherent excitation numbers of the 
photon (a) and the atom (b) fields are plotted as a function 
of the coupling normalized to its critical value yc- With this 
scaling, the mean field amplitudes (thick solid lines) coincide 
in the two states. There is a marked difference between the 
excitation numbers in the steady state (thin solid lines) and in 
the ground state (thin dashed lines). The points show the nu- 
merical results for finite atom number, which we obtained us- 
ing a general quantum simulation framework for cavity QED 
systems |30) . The steady state (circles) and the ground state 
(squares) are calculated for A'^ = 200 and 400, respectively. 
The parameters are fc = — 2, k = 2 con. 
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respectively. Fig. [2] presents a comparison of the quan- 
tum critical behavior in these two cases. The operators 
are split into their steady-state expectation values and 
quantum fluctuations. 



a = VNao + Sa; Ck = e-'^^VN-fk + Sck) 



(4) 



where ao corresponds to a coherent state in the cav- 
ity, and 7o, 71 are the condensate wave function com- 
ponents. The trivial part of the time evolution driven 
by the chemical potential fi is separated from the atomic 
operators. We can choose 70 = y^l — /3g and 71 = (3o, 



FIG. 3: Critical exponents for the phase transition in ground 
state and in the steady state. The scaling of the condensate 
depletion with the relative deviation from the critical point is 
shown in a log- log scale on both sides of the critical point. At 



the center 1 1 — y/yc 



We compare the averages in the 



steady state (solid red lines) for k = 2ijjr and in the ground 
state of the system (dashed blue lines) without dissipation 
(k = 0). The corresponding exponents are —1.00 and —0.50, 
respectively. Other parametes are 5c = —2. 



is independent of k and is related to the different physical 
origins of the divergence. On the one hand, the ground 
state is a two-mode squeezed state with a squeezing pa- 




shape the research on quantum phase transitions, too. In 
particular, we beheve that the cavity QED based systems 
given in several laboratories worldwide raise relevant new 
aspects for classifying quantum critical phenomena. In 
this paper we revealed non-equilibrium critical effects in 
the case of the simplest possible environment. It consists 
of a single well-defined dissipation channel (photon leak- 
age out of the cavity) , which is equivalent with the back- 
action of a weak quantum measurement on the photon 
number observable (continuous photo-detection). This 
innocent looking intrusion in the system drastically mod- 
ifies the critical exponent of the singularity at the phase 
transition point. 



FIG. 4: Entanglement at the critical point. Logarithmic neg- 
ativity as a function of the transverse pumping strength y for 
the steady-state (soUd red lines, k = 2 ljr) and for the ground 
state (dashed blue lines, k, — 0). Parameters: Sc ~ —2. 



rameter which tends to infinity at the critical point. This 
results in a singularity of the entanglement between the 
cavity and the atomic subsystems. On the other hand, 
the steady state is driven by quantum noise associated 
with dissipation (or measurement), which heats up the 
quasi-normal mode population infinitely where the imag- 
inary part of its eigenvalue vanishes. The steady-state is 
a mixed state having a regular entanglement at the crit- 
ical point, reflected by the logarithmic negativity E_\f in 

Fig.m 

The steady-state divergence occurs, as mentioned 
above, when one of the quasi-normal modes has zero 
damping. This means a critical slowing down in reaching 
the steady-state. Therefore, in an experiment where the 
system is typically launched from the quasi-adiabatically 
tuned ground state, the photon number generated in the 
cavity during a finite time is expected in between that of 
the ground state and the steady state. Below threshold, 
where the mean field vanishes, any detected photon cor- 
responds to the fiuctuations (Sa^Sa), and a measurement 
above the ground state level indicates photons generated 
by the quantum noise penetrating into the cavity. 



IV. DISCUSSION 

In this paper we have adapted the famous Dicke-model 
to an intrinsically non- equilibrium setting and pointed 
out distinctive features of this experimentally accessi- 
ble driven-damped open system. The underlying Dicke 
problem, a closed, conservative system, has been a sub- 
ject of intensive research for many decades. It remained 
an intriguing question, however, what happens with the 
critical point under non-equilibrium effects? While in 
classical physics the extension from equilibrium to non- 
equilibrium systems has been extensively studied, this 
step has not been made in quantum theory. On the other 
hand, the ongoing experimental work will significantly 
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Appendix A: Mean-field solution of the two-mode 
model 

We solve the steady state of the system within a mean- 
field approach, which is consistent with the assumptions 
that (i) there is a macroscopically populated BEG wave 
function, and (ii) the state of the cavity field is close to a 
coherent state. After restricting the atomic dynamics to 
the spatial modes Eq. (2) (see also Ref . ^ ) , we proceed 
by obtaining the equations of motion of the operators a, 
Co, ci. The first equation is the Heisenberg-Langevin 
equation given by Eq. (3), while the Heisenberg equa- 
tions of motion for the operators c^ follow directly from 
the Hamiltonian dynamics provided by Eq. (1). 



c[ci 



N 



y 



Co = 1 



C\ 



2 



u 

2N 



a* a 



Co 



y 

2y/N 



(cjci 



cjco 



a' — a) Ci , 



(Ala) 
(Alb) 



2A^ 



a' a 



Cl 



y 



+ ^=(a^-a)co, (Ale) 



2VN 



where we introduce the parameters Sc — Ap — 2u, uji = 
hP/{2m), u = NUo/4: and y = V2N7]t. Note that in the 
thermodynamic limit N ^ oo and Ua,'qt — >■ 0, while u 
and y are kept constant, and they can be expressed with 
the atom density. The noise operator ^ has zero mean 
and its only non- vanishing correlation is {£,{t)^'^{t')) = 
2KS{t — t'), with 2k being the photon loss rate [3ni32| . 

By using the decomposition of the operators to mean 
part and fiuctuations, given by Eq. (4), in the equa- 
tions of motion ( Al I and neglecting fluctuations we arrive 



to the mean- field equations determining Oq, 70 and 71. 
Since the BEC wave function is normalized to unity, we 
can choose 70 — \/l ~ /3q and 71 = /3oi where < /3o < 1- 
By choosing /3o positive, we select one of the two mean 
field solutions, that also fixes the phase of the cavity field. 
For the steady state, we obtain 



[i ((Sc - u^o) - k] "0 = -y/3oY 1 - /^o 



(wfl + u|aoP) /3o = -y Ini(ao) 



1 - 2/?g 



(A2a) 



(A2b) 



and the chemical potential /i = — ^(cjj^ + u|q;oP)/(1 — 
2/30^). 

Note that the coherent field amplitude ao is complex, 
while /3o is real. The solution ag = /3o = always satisfies 
these equations, and it corresponds to the normal phase 
in which the condensate is homogeneous (70 = 1) and 
there is no photon inside the cavity. Above the pumping 
threshold, y\ — —ujji{5'Jj-\-k?)/5c, the solution bifur- 
cates, and the normal phase looses stability. The stable 
solution becomes 



Pi- 



1- Wi- 



u y^ 



Vc 



6c y^ + uujR 



(A3) 



which corresponds to the superradiant phase, where the 
condensate is modulated (/3o > 0) and the cavity fie ld is 
finite (|aoP > 0). For u = 0, the expression in Eq. (A3) 



needs to be reformulated as /3q = {y — y^)/2y^, and 
accordingly the mean field amplitudes coincide both for 
the steady state and for the ground state {k = 0), if 
expressed as a function of y/yc, leaving k the only role 
of shifting the critical pumping strength j/c- 

The critical behaviour is unaffected by u, and the pa- 
rameters u and y can be tuned independently, thus for 
simplicity we set u = in the discussion of this Letter. 



Appendix B: Fluctuations in the steady state 

To go beyond mean field one has to keep the oper- 
ator valued fluctuations Sa and 6ci in Eqs 
consider quantum fluctuations up to linear ore 



(|Alj). We 

Note, 

that the zeroth order term vanishes due to the mean fleld 
equations, and we arrive to a set of linear, stochastic dif- 
ferential equations for the fluctuations. There are two 
types of fluctuations in the atom field (Scg, Sci). The 
zero mode fiuctuations, Sc — \J\ — /SqJcq -I- /3o(5ci, give 
rise to a phase diffusion of the condensate. The dynam- 
ics of the zero mode decouples from that of the other 
types of fluctuations. We are intereste d in the dynamics 
of the non-zero mode 5h = — /3o<Jco -I- \/\ — /3q(5ci, which 
describe the condensate depletion SN — {Sb^Sb). The 



coupled equations of motion read 
d 



dt 



Sa — [i ((5(7 — W/3q) — k] (5a -|- ^ 



(1 ~ 213^^) ~ luaoPo^/^'^ 



[Sb'^ + 5b) , (Bla) 



^B. 



1-2/32 2^ 



2fil) {5a^ - 5a) 
m/3o\/l-/3^(ao'5a^+«oH, (Bib) 



dt 



We solve Eqs. (TBlk-b) by calculating the normal mode 
excitations of the system. Arranging the fluctuations in 
the vector R = [5a, 5a\ 5b, 5b^, Eqs. (Bl i-b) are written 
in the compact form 



d_ 
dt 



R = MR + C : 



(B2) 



where M is the linear stability matrix of the mean field 
solution, and the driving term ^ = [^,^', 0,0] includes 
the quantum noise of the cavity field. The matrix M 
is non-normal, therefore it has differrent left and right 
eigenvectors V- ' and A^\ that form a biorthogonal sys- 
tem, i.e. their scalar product {V- ',7:*^'') = (5^. /. The quasi- 
normal modes defined by pk = {I , R) a-re decoupled 
from each other, and evolve as 



Jo 



Qkit')dt'. (B3) 



Generally, the noise enters into all quasi-normal modes 
via the projection Q^. ~ (r- \^). Since R contains the 
fluctuation operators twice (the operators and their Her- 
mitian adjoint), the operators pk also form adjoint pairs 
P-^- , p\_ with eigenvalues A+ , X*^ and p_ , p_ with eigen- 
values A_,A1. Each pair corresponds to a quasinormal 
mode excitation of the system. Fig. [5] shows the spec- 
trum of the linear stability matrix M as a function of the 
pumping strength y. Solid lines correspond to the imag- 
inary parts, dashed lines to the real parts of the complex 
conjugate pairs of eigenvalues A^, A^. (The real parts are 
the same for each pair.) The conservative BEC and the 
lossy cavity modes form two dissipative quasi-normal ex- 
citation modes, both are subjected to quantum noise and 
damping. One of them is a photon-like mode, with eigen- 
values A+, A^ starting from —k ± iSc at y — 0, with a 
slightly increasing frequency and decreasing decay rate 
as y — >■ yc. The other quasi-normal mode dominantly 
corresponds to the BEC mode. Its eigenvalues A_, \*_ 
are purly imaginary, dziujR at y = 0, however, with in- 
creasing y its decay rate increases and its frequency de- 
creases down to zero. Interestingly, there is a finite in- 
terval where the imaginary part of A_ vanishes. At the 
lower and upper limits of this interval, the matrix M be- 
comes defective, i.e. it has only three independent eigen- 
vectors and three eigenvalues with A_ = A^ becoming a 




FIG. 5: Spectrum of the linear stability matrix M vs the 
transverse pump strength y. Solid lines (dashed lines) correa- 
spond to the imaginary (real) part of the eigenvalues. The real 
parts are the same for a complex conjugate pair. Parameters: 
Sc = —2, K = 2 ujR. 



multiple eigenvalue. The critical point is reached, where 
the smallest decay rate becomes zero. At this point, the 
quantum noise is not balanced by damping, therefore the 
steady-state excitation numbers diverge. 

The second order correlations of the original fluctu- 
ation operators can be derived from the correlations 
{Pk{t)pi{t)). In the regime of cavity cooling, where 
Sc — uf^Q < 0, the real parts of the eige nvalues Aj, are 
negative, thus the first term of Eq. (B3| dies out with 
time. The steady-state correlations are then obtained 
from the second term 



{Pk{t)pi{t)) 



2k 



,(fc)*,(0* 



^k + A/ 



(B4) 



This result is in contrast to the normal mode expecta- 
tion values of zero temperature systems. For such an 
equilibrium situation the expectation values are simply 
(pkPi) = 1, provided that pk is the annihilation, and pi 
is the creation operator of the same normal mode, i.e. 
< Im(Az) = - Im(Afc). 

The correlations of the original system operators can 
be calculated using their expansion with the quasi normal 
modes, R — ^^ Pkn'^'^K that leads to 



(i?4j) = Y.(p'^pi) 

k,l 



I 3 



(B5) 



For expample, the condensate depletion is given by 5N = 
{db^db) = {R4R3), while the number of incoherent cavity 
photons is expressed by (Sa'^Sa) ~ {R2Ri). 



Appendix C: Atom-field entanglement 

We quantify the entanglement between the EEC and 
cavity subsystems by calculating the logarithmic nega- 
tivity from the steady-state correlation matrix. To this 
end, we introduce the quadrature operators Sx = {Sa + 
(5at)/V2, Sy = -i{Sa - 5a))lV2, 5X = {Sb + Sb'')/V2, 
6Y = —i{Sb — db^)/y/2 and group them in the vector 
u = {Sx,6y,6X,6Y)'^ . As these quadratures are hermi- 
tian, one can construct a real correlation matrix 



Cij = ^{uiUj +UjUi), 



which has the following block form 



(CI) 



c = 



P X 

X^ A 



(C2) 



where P and A describes the correlations whithin the 
photon and atom fields, while X accounts for the cross 
correlations between the two. The logarithmic negativity 
can be expressed |33j by the symplectic invariants (det P, 
det A, det X) of the covariance matrix (CI ) as 



Ej^ = max(0, — log2i>_). 



(C3) 



and 



I](C) - VS(C)2-4detC, 



(C4) 



where S — det P + det A — 2 det X. The state is separable 
thus the entanglement is zero if :>_ > i. The logarithmic 
negativity quantifies the amount by which this separabil- 
ity criterion is violated. 
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